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Abstract 

Many models for physics beyond the Standard Model predict lepton-flavour violating 
decays of charged leptons at a level which may become observable very soon. In the 
present paper we investigate the decays of a r into three charged leptons in a generic 
way, based on effective-field-theory methods, where the relevant operators are classified 
according to their chirality structure. We work out the decay distributions and discuss 
phenomenological implications. 



1 Introduction 



Lepton flavour violation (LFV) has become a hot topic over the last few years. On the one 
hand, the discovery of neutrino oscillations in combination with the minimal extension of the 
Standard Model (SM) predicts lepton flavour violation for the charged leptons, however, at a 
completely unobscrvable level. On the other hand, many extensions of the SM predict LFV at 
much higher rates, which, in some cases, may already be in conflict with existing experimental 
bounds [1] (see also [2] for a recent summary of 5-factory results). With the advent of new 
experimental facilities [3] (see also [4]) the current bounds will be pushed further, if not a 
discovery will be made. In particular, at the LHC experiments it will be possible to detect LFV 
decays of a r lepton, especially into channels with three leptons; here the signal r — > 3// will 
be one of the cleanest signatures [5]. 

There are many models which predict LFV r decays of the form r — > with £, ^" = e, n 
[6-17]. All these models will eventually match onto a set of local four-fermion operators or 
radiative operators, the latter mediating r — > £7* with subsequent decay of the (virtual) photon 
into a charged lepton pair. In a bottom-up approach, this allows us to consider all possible four- 
fermion and radiative operators with arbitrary coupling constants, which can be determined by 
studying the decay distributions of the three leptons in the final state. Even if no signal events 
are found, such a study of the decay distributions is necessary to determine the efficiency of an 
experiment and hence to extract reliable limits. 

Our paper is organized as follows. In the next section we will introduce the general set 
of effective operators of dimension six and eight, which mediate the decays r ii'i" at the 
electroweak scale. This includes four-fermion operators with scalar, vector and tensor currents 
as well as radiative operators contributing to r — > £7 ^ i£'£". These operators match onto the 
relevant low-energy interactions at the scale of the r mass, which are parameterized by a num- 
ber of unknown coupling constants. We will also give a brief discussion on how the couplings are 
related to the lepton masses and the PMNS neutrino mixing matrix within minimal-flavour vio- 
lating scenarios [18]. Focusing on the decays r~ — > ii^/i^ii'^ in section 3.1 and r~ — > e~ ii~ 11'^ in 
section 3.2, we calculate the Dalitz distributions for the individual chirality structures appear- 
ing in the effective Hamiltonian, taking into account interference terms apart from corrections 
of order m^/m,-. We conclude with a brief comparison of existing results on LFV processes in 
specific new physics scenarios in section 4. 

2 The Effective Interaction for r Wl" Decays 

In this paper we follow a bottom-up approach to lepton-fiavour violating decays, using an 
effective-field-theory picture. We assume that some new physics at a high scale A induces 
lepton-flavour violating processes.^ At the electro-weak scale these lepton-flavour-violating 
interactions manifest themselves in higher dimensional operators which have to be compatible 
with the SU{2)l x U{1)y gauge symmetry of the SM. The leading operators for processes 
involving charged leptons will be of dim=6 and dim=8, which we shall construct in the following 
(similar analyses within the context of supersymmetric extensions of the SM can be found, for 
instance, in [8, 10]). 

To this end we group the left-handed leptons in an SU{2)l doublet, while the right-handed 
charged leptons (which are singlets under SU{2)l) are put into an incomplete doublet, as a 

^Notice that, in general, the scale associated to lepton-flavour violation is independent of the scale related 
to lepton- number violation, Aln- 



1 



reminiscent of a right-handed SU{2)ji related to custodial symmetry. Writing also the Higgs 
boson in matrix form, we have 

L=f^^V i?=fM, H^^i'^^^'^J^' V (2.1) 

V y ) a/2 V -^(t>- v + ho-ixo J ^ ^ 

For simplicity, we have suppressed the family indices, which will be specified once we consider 
a particular decay mode. For later use, we note that the hypercharge Y is defined in terms of 
the right-handed generator T^^^ as 

y^rf)+l(S-L). (2.2) 
In terms of the fields defined in (2.1), the list of operators relevant for r — > ££'£" decays reads 



dim=6 leptonic: 

01 = {Lj.LKLYL) (2.3) dim=8 leptonic: 

02 = CLr^l,L)CLr^YL) (2.4) P, ^ {^LH R){LH R) (2.9) 
0^ = {R^^R)iRYR) (2.5) P^^{Lt''HR){Lt''HR) (2.10) 
0, = {R^,R){LrL) (2.6) Q,^{LHR){RH^L) (2.11) 

Q2^{Lt''HR){RH'^t''L) (2.12) 

dim=6 radiative: 

R,^g'iLHa„.R)B^ (2.7) P^^' = (LHa^mLHa'^R) (2.13) 

ft, = 9(Lr°i/a^fi)W"«'.» (2.8) PP ^ (Lt'H^^R)(It'H,jI-R) (2.14) 



where B^^^ is the U{1)y gauge field, W^^ are the SU{2)l gauge fields, and g and g' are the 
corresponding gauge couplings. In the above list wc have only shown operators that have tree- 
level contributions to leptonic r-decays; more operators, which are bi-linear in the lepton fields 
and contribute at the loop level, can be found e.g. in [8, 10, 18]. We also neglected dim=8 
operators involving additional covariant derivatives. When acting on fermions, the derivatives 
become fermion masses by the equations of motions, such that these operators are additionally 
suppressed by the small lepton Yukawa couplings. Notice that the leptonic dim=6 operators 
O1 — O4 only contain helicity conserving currents. The most general effective Hamiltonian at the 
electro-weak scale is then obtained by summing over these operators, multiplied by arbitrary 
coefficients for every fiavour combination. In a particular new physics scenario, these coefficient 
should be obtained by matching at the new physics scale A and evolving down to the scale Mw 
within the SM (as an effective theory). 

In the following we are interested in LFV decays of a r lepton into three charged leptons. To 
this end we have to construct the effective interaction at the scale of the r lepton, by integrating 
out the weak gauge bosons and the Higgs. We will focus on r" decays; the decay distributions 
for r"*" decays are identical. On the level of four-fermion operators with dim=6, we obtain the 
same structures as in (2.3-2.6). Projecting on charged leptons only, we see that O2 becomes 
equivalent to Oi, and both match onto a purely left-handed operator 

rriLL)iLL) _ (LL){LL) {^Ll hTlWlY ^-'D . 
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where here and in what follows the superscript of the coupling denotes the combinations of chi- 
ralitics involved and the subscript denotes the relevant Dirac structure. Likewise, the operator 
O3 corresponds to a purely right-handed interaction 

rr{RR){RR) _ (RR){RR) {^RIhTr) {^'Rl^f-R) , . 

-"eff — 9v 'J^ ' \'^-^^) 

while we get a mixed term from the operator O4 

rr{LL){RR) _ (LL)(RR) {^11 ^,Tl){1'rY ^r) , (RR){LL) {^RlinTR){i'^'y 

Notice that the chirality structure (LR) (RL) is not independent, since it can be Fierz rearranged 
intoi/(i^)(^^). 

The dim=6 radiative operators contain charged as well as neutral currents. Here we are 
only interested in the neutral-current component, coupling to a charged lepton pair. Switching 
to the physical photon and Zq fields and integrating out the Higgs and the Zq, wc find that 
we obtain a radiative operator with a photon, as well as a four-fermion contribution from Zq 
exchange. The latter is proportional to 

where (gA) gv are the (axial)vector couplings of the Zq to the leptons. Taking into account the 
fact that Ig'*! is of the order of the r mass, we find that this operator is suppressed relative 
to the leading ones by the small Yukawa coupling of the r lepton. Thus, only the photonic 
contribution has to be taken into account. For this we obtain at the scale itIt 

^l^-JZT2ll ^^^'^ (4(-^<T,.)r.) , (2.18) 



h.s 



where g^^^^ and g^^^^ denote the two possible chirality combinations.^ The matrix element 
for T — > £i'i' becomes 



Hlt\r) =c^em^^Yl ^rt'^ (^^"^'1 {iYi') \r) , (2.19) 



h.s 



where q is the momentum transfer through the photon. This momentum transfer is proportional 
to the lepton masses, and thus this contribution scales as l/(yA^) where y is a Yukawa coupling 
of the leptons, which would lead to an enhancement unless an additional Yukawa coupling 
appears in the numerator as e.g. in minimal flavour violation (see below). 

Finally, we turn to the four-fermion operators with the chirality structure {RL){RL) or 
{LR){LR). At tree-level, they receive contributions from the dim=8 operators Pi,2, Qi,2 and 



^Here, for simplicity, we neglected possible form factor effects for decays into virtual photons from long- 
distance lepton or quark loops. In the most general case, the r ij* vertex could be parametrized as 



e V 
4^ A2 



h,s ^ \ / J 

where gl^J'\o) = ^rad*^ and f^^\o) = 0, see, for instance, [19]. 
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Pi^ , only. Therefore, their matching coefficients are further suppressed by v^/A^. The one- 
loop matching coefficients at the electroweak scale may also receive contributions from the 
dim=6 operators Oi_4, i?i,2, but in this case the required chirality flips induce an additional 
suppression by mj/v^. Ignoring contributions suppressed by v^/A"^ or m^/f^, this reduces the 
number of possible Dirac structures already to six in the case where the radiative operator can 
contribute and to four in cases like t" — > ii~ii~e'^, where the radiative contribution is absent. 
The corresponding couplings 

{LL){LL) (RR){RR) (LL){RR) (RR){LL) (LR) (RL) 

yV ) yy ) yy ^ (Jv ' i/rad ' i/racl ' 

are matrices in lepton flavour space. There are in total six different decay modes of the r~ to 
consider, 



r" ^ e-e-e+ (2.20) 1^-^-6+ (2.23) 

T" i^-^'i^^ (2.21) r" /i"e"e+ (2.24) 

r- e'e-fi^ (2.22) r' e' (2.25) 

Notice that (2.20 - 2.23) contain two identical particles (e~e~ or /i~/i~) in the final state, 
whereas (2.24 + 2.25) do not. Moreover, only (2.20, 2.21, 2.24, 2.25) receive contributions from 
the radiative operators (2.18) via 

2.1 Constraints from Minimal Flavour Violation 

The fiavour structure of the coupling constant has been investigated within the framework of 
Minimal Flavour Violation in the lepton sector (MLFV [18]). In a scenario with the minimal 
field content (2.1), the breaking of the lepton fiavour symmetry SU(i)L x SU{?>)eii described 
by two spurion fields 

A = — = -diag(me,mu,m^) , (2.26) 

V V 

Qv = ^^t/*diag(m^,i,m^2,m^2)t/^ (2.27) 



V 



where A ~ (3, 3) describes the SM Yukawa couphngs of the charged leptons, and the matrix 
Qi, ~ (6, 1) stems from a dim-5 lepton-number violating term. 



-Maj 



2Aln 
where 



{N'^gN) , (2.28) 



N = yr^^^ + H^L (2.29) 

has vanishing quantum numbers under the complete SM gauge group. If the scale Aln, 
associated with lepton-number violation, is sufficiently large, the resulting neutrino masses 
?7iMaj ~ ^^^/Aln are small, even if the spurion g,j has generic entries of order unity. 

We are interested in 4-lepton processes, induced by operators with some fiavour structure 

U U LI L*i , U LI Rl , etc. 
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To render these operators formally invariant under the flavour group, they have to be multiplied 
by appropriate factors of Ag and Qi,. In the following we will consider the minimal number of 
spurion insertions, only.^ 

Starting with the case of U U L^L^, we need at least two spurion insertions. The possible 
flavour structures can be read off the reduction of the SU{3)l tensor product for g^, and gl, 

6x6 = 1 + 8 + 27. 

Here the flavour singlet term corresponds to the trace of glg^, 

^^[9l9A - ^ {< + < + O ^ 3 , 
which does not induce flavour transitions at all. The flavour octet term is obtained as^ 

A = At = glg^-U,[glgJ^ = ^U^mlUK (2.30) 

Here Am^ = diagfm^^, m^^; "^^3] ~ '^t- particular, one flnds that 



^4 ^""atm 

whereas A^ and A^ are further suppressed by the neutrino mixing angle ^13. It is to be stressed 
that A docs neither depend on the absolute neutrino mass scale m^, nor on potential Majorana 
phases ai^2 in the PMNS matrix. 

The flavour structure of the corresponding invariant 4-lepton operator reads (L* A L) {L* L) . 
The coefficients of specific fiavour transitions are thus given by the quadratic neutrino mass 
differences and PMNS elements. By the same argument, the operator (L* A L){R* R) is invari- 
ant under the flavour group. It has been shown in [18] that A also drives all possible flavour 
transitions induced by operators that are bilinear in the lepton fields. In particular, the fiavour 
structure of the radiative operators in (2.7,2.8) reads (L* AA^i?) and (i?*A AL) . Notice that 
the presence of a single right-handed field requires the insertion of the Yukawa spurion A, which 
leads to an additional suppression factor m^/v. 

Turning to the 27plet combination of g^ and g]^, we introduce the according representation 
in terms of a trace-less tensor 

= {g.)^A9lr~Y2^8^ + ^,5])i.{g^g) 

- \ {5%5^, + ^'A^i + ^t^A + ^A^^ K (2-31) 

with Gfj = Gjl = Gl'j, and Gfj = 0. The fiavour structure of the corresponding invariant 
4-lepton operator reads 

-ikl 



^This is justified as long as the spurion fields are characterized by some small expansion parameter [18,20,21], 
e.g. if the neutrino mass differences Ami smaller than their average Affif,. Notice that, unlike in the case 
of the quark CKM matrix, the off-diagonal entries of the PMNS matrix are not always small. 

^Notice that our definition of A differs from the one in [18], but only for the diagonal elements, which are 
irrelevant for flavour transitions. 
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In contrast to A, the off-diagonal matrix elements of G depend on the absolute neutrino mass 
scale fhl and the Majorana phases. As a consequence, in the general case, i.e. if the radiative 
operators do not dominate the r ^ 3£ decay amplitudes, the purely leptonic decay modes 
are not directly correlated with the radiative ones r — > £7. Relatively simple expressions for 
Gfj can be obtained in the limit of vanishing Majorana phases, where we also employ the 
approximations sin^ 6'i3 ~ Amg^^/Am^^j^ -C 1 and 623 = 45°. For the normal neutrino hierarchy 
{rrii.^ ~ m,^2 <S "^i^s), we obtain the leading coefficients as 

GTe ^ -2G% , (2.32) 

and sub-leading effects from 



— e sm Am^tm , ^r/,--^ ^ ( ^^^1,2 - cos 2^12 — — 



1^1,2 



C.e - [e sm e,, V Am,,^ - sm 29,, ^ J , 

C., ^ ^ ^ ^ sm VAm,,^ + sm 2^1^ ^ J • (2.33) 

where we have used the PDG parameterization [1] of the CKM matrix. For the inverted 
hierarchy {niy^ ~ rUy^ 3> rriy^), one has 



Gf^ ~ -5G^^ ~ lOG^!^ ~ ^ — ^ , (2.34) 



and 



Aln Ami^ . . 3 cos (5 - 7i sin (5 

G^e - ^ sm ^,3 ^ , 

~ sin ^13 , Gt^^^^ ^^^^ e^^ sin ^13 ■ (2.35) 

We finally note that purely right-handed lepton-flavour violating decays require at least four 
spurion insertions, {R*R){R*Xg'^g X^R) , and are thus strongly suppressed in MLFV. 

In summary, to obtain the dominating flavour coefficients for the operators in (2.15, 2.17, 
2.18), relevant for flavour- violating r decays in MLFV, one has to consider 

giL.L^KULn ^ 2ciAf5j. + C2Gg, (2.36) 

giL.L^mR^) ^ ^^^,^1^ (2.37) 

dlT^ - C4Af, (2.38) 

whereas the chiral structures corresponding to g^^^^^^^ ^ g{RR){LL) gl^^"* are suppressed by 
small lepton masses. The spurion combination G^j represents a new source of LFV compared 
to the radiative transitions r — > ^7. While the latter are driven by the spurion A and hence 
by the difference of the squared neutrino masses, the flavour coefficients of purely left-handed 
four-lepton operators in MLFV also involve the absolute neutrino mass scale as well as the 
Majorana phases. In particular, the decay modes (2.22 -|- 2.23) only receive contributions from 
and G:^^^, respectively. 
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Figure 1: Dalitz plot for d^r^^^^^^^^ (left) and dT^; 
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3 Dalitz-Plot Analysis 



3.1 The Decay r 



In this section we will give a detailed analysis of t~ — > ii~ ii" ii^- as the probably most prominent 
channel to be looked for at the LHC To this end, we shall consider the Dalitz distributions for 
the different chirality structures (2.15, 2.16, 2.17, 2.18) in the dim=6 effective Hamiltonian, in 
terms of the variables 



m 



m 



12 



= (P^-+P;-)', 



m 



+- 



m 



23 



(3.1) 



and mfg = + 3m^ — — m^_. We will make use of (approximate) helicity conservation, 
which implies that many of the interference terms between the operators with different chiralities 

are suppressed by powers of m^, and can be ignored to first approximation. 

In the simplest case all four leptons are left-handed, and the decay amplitude is determined 
by H^Q^^''^^^ in (2.15). The corresponding Dalitz distribution 



2p(LL)(LL) 



d^r 



V 



dm^g dm^2 



I (L^L-)(L^LM) |2 

19V \_ 



{ml — m?^Y — (2m 



^12 



2\2 



3m^) 



256 TT^ 



(3.2) 



is shown in Fig. 1 (left). The events are equally distributed along m^_, while there is a rather 
flat maximum at 



m 



12 



analogous and yields the same distribution with (?[^^^^^^^ 



ml/2. The case with all particles right-handed is completely 

^{jriR)(RR) ^ (We remind the reader 
that g^^''^^^'' is expected to be strongly suppressed within MLFV scenarios.) 

Next we will consider the operator ff^^^^^^^^ in (2.17). For a left-handed r-lepton we obtain 
the Dalitz distribution 



2j.(LL){RR) 



d"r 



V 



dm23dmf2 



|^(L^I,-)(Rm-R'')|2 
A4 



{ml — m^^Y — 4m^ {ml + — m\2) 



512 TT-^ 

2 ^2 



.2\2n 



(2m^3 — ml — Smj^y + (2m23 — ml — 3mjJ 
1024 n^ml 



(3.3) 



shown in Fig. 1 (right). In this case the events are distributed around a flat maximum at 
~ ml/2 and ~ 0. Again, the case of a right-handed r yields the same distribution. 
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As pointed out above, the interference terms between (2.15) and (2.17) are suppressed by 




m\_ I GeV^ 

Figure 2: Dalitz plot for d^F^^^-* in t~ — > ix~ ix~ [x^ . 



In addition to the four-fermion operators we also get a contribution from the radiative 
interaction via (2.19). The resulting Dalitz distribution for a right-handed r-lepton, 



d"F 



rad 



dmiadmfa 



\9- 



(iM-R")|2„,2 



rad 



V 



m 



I {ml - mlf 
1287r3m? 



1 1 

+ 



rri 



13 



23 



+ 



128 TT^ mfg 77123 '^r 



(mfa + mi3)(mf2 + mf^ + mj^ - 6m^(m^ + ml)) 2m?2 - 3m| 

~t~ ^ ^ O O O Q ~t~ 



256 TT^ 77123 "^r 



1287r3m3 



(3.4) 



is plotted in Fig. 2. Due to the photon pole, the events are concentrated at low values of 7?7|3 
or 771^3, respectively. Again the decay of the left-handed r is completely analogous. 

Finally, we have to take into account the contributions from the interference terms between 
the radiative operators and four-fermion operators, for the cases where only the chirality of the 
r-lepton has to be flipped. The interference term between (2.15) and (2.18) reads 



d^F 



2-^(.LL){LL) 
mix 

dml^dml2 



2 V tie[gy g^^ J 



A4 



64 TT^ 771: 



+ 



128 TT^ 777^ r?7i3 m^^ 



The interference between (2.17) and (2.18) results in 



(3.5) 



dT 



2-^(.LL){RR) 
mix 

dr?7i3dr?7f2 



A4 

ml - ml^ - Sml ml(ml - mDimj^ + mig) 

X :^ r 



256 TT^ r?7^ 



256 TT^ r?7^ r?7f3 m^^ 



(3.6) 



In both cases, the photon pole at 777^3 = or 777I3 = is suppressed by the small muon mass. 
The remaining terms increase (decrease) monotonically with 777^2, respectively, see Fig. 3. 



8 



m?._ / GeV^ ml_ / GeV^ 

Figure 3: Dalitz plot for \<l^T^ii!^^^^\ (left) and Id^r;^^^^^^^^ | (right) in ^ 11-1^-11+. 

Combining (3.2-3.6) and integrating over phase space, we obtain for the total decay width 
(normalized to the SM decay r — > /iP^iyr, and neglecting the muon mass) 

r[r- ^ ^ 1 / 2 + 2 l^r^^''''^ + [^{.'^^^^'^'^^ + 

r[T- ^ /.-i/^i/,] G|A4 I 8 

Oen^ r *(LR) (LL){LL) *{LR) {LL){RR) , / r ^ px" 

This result is consistent with the formula quoted, for instance, in [10]. 
3.2 The Decay r" e' ii' 

For completeness, we also discuss the Dalitz distributions for the decay mode e~ij,~ij,'^, 
which belongs to the class of decays with three different particles in the final state. Again, we 
will give the results in terms of the invariant masses 

m^__ = ml^ = (Pe- +P'„-)^ , rnl_ = ml^ = {p'^- +p^+f , (3.8) 

and = + 2m^ — — m^_, where we set the electron mass to zero.^ 

From the purely left-handed term in the effective Hamiltonian, H^q^''^^^^ in (2.15), we obtain 
the Dalitz distribution 

which (except in the vicinity of the phase-space boundaries) coincides with (3.2) up to cor- 
rections of order m^/m^ and a statistical factor. Consequently, the corresponding Dalitz plot 
looks almost identical to Fig. 1 (left). 



^Notice that the photon pole from r — > e 7* — > e /i /U+ is still regulated by the muon mass. 
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m?._ / GeV^ 



ml_ I GeV^ 



Figure 4: Dalitz plots for the two contributions to d 



2-r^(LL){RR) 



m T 



Prom in (2.17) we obtain for the case of a left-handed r-lepton 

(2m^3 - - 2m^)^ 



^2p(LL)(Ril) 



\9v l_ '"'T 

(L^L^)(i?eR'')|2 



+ 



\9'v 



m; - 
512 TT^ ml 

,2 o^2\2 



2\2 



,m; — 2m^)^ — (2m23 — f^l — 2rrr^) 



512 TT^ m2 



(3.10) 



The corresponding Dalitz plots for the two separate contributions are shown in Fig. 4 (left: the 
term oc {g^^"^ )|2. j^jgj-^^- ^]-^g term oc \g\^^^ )iReR )^2y rjj^g events are distributed around 

+ ^ m^/2 or ~ fn^/'^-, respectively. For equal coupling constants in (3.10) we 
recover the r — * 3/i case in (3.3) (again up to mass corrections and a statistical factor). 
For the radiative decay operators, we obtain 



j2p(Li?) 
^ rad 

dmlgdm^s 



If, 



(iei?-)|2^2 



rad 



ml (m^3 - mlf 

64 TT^ 777.23 



12 



13 



2777^ 



128 TT^ 777^ 777^3 



+ 



mZ — m 



23 



128 7r3 7773 



(3.11) 



and the corresponding Dalitz plot is shown in Fig: 5. In this case the photon pole enhances the 

events at low values of 7r7i 



m: 



23- 



Finally, for the interference terms between (2.15) and (2.18) we get 



.2j.{LL){LL) 
^ mix 



d77723d777? 



= Clip 



2 7;Re[,{f^-)(^''^''),;i^-)] 



12 



A4 



777 



12 



2777^ 



1287r3777? 



+ 



7777 



1287r3777 ^3 



(3.12) 



and 



2-piLL)iRR) 



dml^dmj^ 



2 7;Re[5{f^^^)(^''^^)^;i^^^^)] 



A4 



7r7 



13 



2777?, 



VlM'^ml 



+ 



m, 



1287r377723 



(3.13) 



The corresponding Dalitz plots are shown in Fig. 6. 
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4 Discussion and Conclusions 



Lepton-number violating processes, like the decay r~ — > discussed in this paper, provide 
an important test of the Standard Model. Since many new physics models allow for dramatic 
enhancements compared to the tiny SM effects, there is a potential to falsify the SM by mea- 
suring such decays at future experiments. At the same time, the foreseen improvement of 
experimental limits may further tighten the constraints on specific new physics models. 

In both cases, the event distributions in phase space provided in this work, will be helpful. 
The difference between the Dalitz distributions arising from four-lepton operators, which show a 
rather uniform behaviour, and the distributions from radiative operators which are concentrated 
at small values of = clearly gives a handle to disentangle different new physics 

models already on the basis of rather few events. Typically, these models give rather different 
predictions for the relative size of radiative and four-fermion operators, as we will discuss in 
the following. 

In super-symmetric extensions of the SM, one typically finds that the photon-dipole oper- 
ator, induced by penguin diagrams, dominates over the four-lepton operators. This leads to 
simple correlations like (see e.g. [10]) 

^^^^ fin 4-^1 =(^(10-^) • 

In this case, one expects Dalitz distributions as shown in Fig. 2. 

It has been pointed out in [13] that Higgs-mediated r — > transitions may alter this result, 
if tan j3 and the off-diagonal slepton mass-matrix element are large. For instance, in the 
decoupling hmit (cos(/3 — a) = 0, m^o » Mz) the author of [13] finds 

r(T ^ ^^x^J) 3 + hSt^, 

r(r - ^7) - ^ ^ ^ 

where = mj^/m^. Testing such scenarios in experiment will be more involved, as one 
generally has to allow for the interplay of all contributions to the Dalitz distributions. Fig. 1-3. 

The situation is somewhat different in the case of Little Higgs Models with T-Parity 
(LHT) [15-17]. Here the Zq and box-diagram contributions dominate compared to the radiative 
operators [17], which is mainly due to the constructive (respectively destructive) interference 
between the individual heavy gauge boson contributions. Depending on the parameter values 
of the LHT, one finds 

3^ = o(i) 

r(T ^ /it) ^ ' 

for a mass scale of the LHT mirror fermions of about 1 TeV. In this case, one can expect rates 
for LFV decays which are already close to the present bounds. Because of the sub-dominance 
of the radiative dipole operator, we expect a rather fiat Dalitz distribution for r — > 3/x, as 
illustrated in Fig. 1. 
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